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five  independent  internal  variables  and  two  connected  feedback  circuits,  either 
positive  and  negative  circuits  of  lenght  4  or  positive  circuits  of  lengths  4  and  2. 
The  flow  topology  of  either  of  these  systems  indicates  that  addition  of  the  fifth 
vertex  switches  the  basic  foui^vertex  core  regulatory  circuit  into  one  of  its  two 
stable  configurations. 
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Immunological  control  networks  can  be  modelled  by  bipartite  signed  directed 
graphs  called  influence  diagrams.  Possible  flow  topologies  around  unstable  steady 
states  in  such  systems  can  be  determined  by  kinetic  logic  based  on  the  directed 
Boolean  cubes  of  switching  circuit  theory.  This  method  can  identify  the  stable 
"help"  and  "suppression"  configurations  of  the  basic  Herzenberg  core  regulatory 
circuit  as  well  as  the  "memory"  and  "non-responsive"  configurations  of  the  model 
of  Kaufman,  Urbain,  and  Thomas.  Furthermore,  addition  of  a  single  vertex  to 
the  Herzenberg  core  regulatory  circuit  so  as  to  preserve  the  bipartite  nature 
of  the  signed  directed  graph  leads  to  a  network  with  five  independent  internal 
variables  and  two  connected  feedback  circuits,  either  positive  and  negative  circuits 
of  length  4  or  positive  circuits  of  lengths  4  and  2.  The  flow  topology  of  either 
of  these  systems  indicates  that  addition  of  the  fifth  vertex  switches  the  basic 
four-vertex  core  regulatory  circuit  into  one  of  its  two  stable  configurations. 
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1.  Introduction 

A  research  area  of  major  interest  to  immunologists  is  the  nature  and  extent 
of  the  interactions  that  regulate  the  quantity  of  antibody  secreting  cells  produced 
in  response  to  antigen.  In  this  connection  network  theory  has  been  used  by  a 
number  of  workers  to  model  such  interactions.  The  general  objective  of  this 
pacer  is  the  examination  of  such  immunological  control  networks  using  topological 
and  graph-theoretical  methods  similar  to  those  previously  used  for  chemical 
reaction  networks  (King,  1980,  1982,  1983).  Such  methods  are  based  on  kinetic 
logic  (Thomas,  1979)  and  in  the  case  of  immunological  control  networks  use  the 
flow  topology  in  Boolean  or  logical  variable  space  to  identify  which  of  the  proposed 
networks  can  exhibit  stable  configurations  corresponding  to  the  experimentally 
observed  "memory"  and  "non-responsive"  states  (Kaufman,  Urbain,  and  Thomas, 
1985).  The  extension  of  these  methods  from  chemical  to  immunological  control 
networks  is  not  trivial  since  the  general  patterns  of  immunological  control 
networks  are  quite  different  from  those  of  chemical  reaction  networks. 
Immunological  control  networks  treated  in  this  paper  include  those  which  can 
be  modelled  using  no  more  than  five  variables.  Such  networks  include  networks 
proposed  by  Herzenberg,  Black,  and  Herzenberg  (1980);  Kaufman,  Urbain,  and 
Thomas  (1985);  Hiernaux  (1977);  and  Richter  (1978a,  b).  Other  networks  such 
as  the  "plus-minus"  network  of  Hoffmann  (1978)  require  a  larger  number  of 
variables  and  interactions  and  therefore  are  not  as  conveniently  treated  by  the 
methods  outlined  in  this  paper. 

Two  research  groups  have  recently  published  work  on  immunological  control 
networks  using  methods  related  to  those  discussed  in  this  paper.  Kaufman,  Urbain, 
and  Thomas  (1985)  have  used  a  similar  but  not  identical  flow  topological  analysis 
also  based  on  the  ideas  of  kinetic  logic  to  evaluate  one  of  the  models  also  discussed 


in  this  paper.  Eisenfield  and  DeLisi  (Eisenfeld  and  DeLisi,  1985;  Eisenfeld,  1986) 
have  analyzed  the  qualitative  stability  of  immunological  control  networks  using 
network  diagrams  to  depict  essential  topological  relationships  and  mathematical 
methods  related  to  those  used  by  Clarke  (1980)  for  the  study  of  chemical  reaction 
networks.  This  paper  extends  the  approach  of  Eisenfeld  and  DeLisi  by  considering 
the  flow  topologies  rather  than  the  stabilities  of  immunological  control  networks 
and  treats  a  greater  variety  of  immunological  control  networks.  Such  topological 
models  provide  guidance  for  much  more  complicated  continuous  analyses  using 
specific  systems  of  differential  equations. 


2.  Background 


An  immunological  control  network  can  be  represented  by  a  bipartite  signed 
directed  graph.  The  vertices  of  this  graph  represent  distinguishable  subsets  of 
lymphocytes.  The  lymphocytes  encountered  in  the  models  in  this  paper  are  helper 
and  suppressor  thymus- derived  (T)  lymphocytes  (Th  and  Ts,  respectively)  and 
the  so-called  B  lymphocytes.  Each  lymphocyte  has  an  associated  idiotype  which 
may  be  positive  or  negative  corresponding  to  a  coloring  of  the  vertices  of  the 
graph  in  one  of  two  colors.  In  the  graphs  used  in  this  paper  vertices  associated 
with  positive  and  negative  idiot-ypes  are  represented  by  squares  or  rectangles 
and  circles  or  ellipses,  respectively. 

The  directed  edges  of  the  immunological  control  network  represent  interactions 
between  lymphocytes  in  different  subsets  based  on  complementary 
idiotype-anti idiotype  interactions  (Jerne,  1974).  Thus  edges  can  be  directed 
only  from  a  positive  to  a  negative  vertex  or  from  a  negative  to  a  positive  vertex. 
Edges  between  pairs  of  vertices  of  the  same  sign  (i.e.,  positive  to  positive  or 
negative  to  negative)  are  not  possible.  An  immunological  control  network  can 
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thus  be  represented  by  a  bipartite  graph.  The  directed  edges  are  given  positive 
or  negative  signs  depending  upon  whether  they  represent  help  or  suppression, 
respectively. 

The  vertices  of  an  immunological  control  network  can  be  classified  as  either 
strong  or  weak  vertices.  A  strong  vertex  has  at  least  one  edge  directed  towards 
it  and  at  least  one  edge  directed  away  from  it,  e.g. 

— KO — * 

I 

A  weak  vertex  either  has  only  edges  directed  towards  it  (a  sink,  lla)  or  only  edges 
directed  away  from  it  (a  source,  lib), 

lib 

Strong  vertices  represent  internal  variables  and  weak  vertices  represent  external 
variables.  A  strong  vertex  with  two  or  more  edges  directed  towards  it  is  called 
a  turbulent  vertex. 

The  relationships  between  the  internal  variables  are  of  interest  in  the  treatment 
of  immunological  control  networks  in  this  paper  as  well  as  in  the  previous  work 
on  chemical  reaction  networks  (King,  1980,  1982,  1983).  Such  relationships  are 
represented  by  signed  directed  graphs  known  as  influence  diagrams.  An  influence 
diagram  is  obtained  from  the  corresponding  immunological  control  network  by 
deleting  all  of  the  weak  vertices  and  any  edges  connected  to  weak  vertices  until 
only  strong  vertices  remain.  Such  deleted  weak  vertices  represent  external 


variables  such  as  antigens  (sources,  e.g.  lib)  or  antibodies  (sinks,  e.g.  I  la). 
Conversion  of  an  immunological  control  network  into  the  corresponding  influence 
diagram  can  require  several  stages  if  removal  of  a  weak  vertex  and  its  associated 
edges  converts  a  strong  vertex  into  a  weak  vertex. 

This  procedure  is  illustrated  by  the  conversion  of  an  immunological  control 
network  reported  by  DeLisi  (1983)  to  the  corresponding  influence  diagram  (Figure 
1).  In  the  original  network  (Figure  1,  top),  the  vertices  A  and  B  are  weak  vertices 
and  SI,  S,  HI,  H  and  M$  are  strong  vertices.  Removal  of  the  vertices  A  and 
B  and  the  associated  four  edges  A  ^  SI,  A+  HI,  H-*-  B,  and  M$-»-  B  leads  to  a 
network  with  the  five  vertices  SI,  S,  HI,  H  and  M$.  (Figure  1,  middle);  all  of 
these  five  vertices  are  strong  vertices  as  required.  In  this  case  the  relationship 
between  the  circuits  in  this  diagram  can  be  represented  by  a  graph  (e.g..  Figure 
1,  bottom)  in  which  both  vertices  are  weak.  This  graph  represents  the  relation 
between  the  suppressor  T-cells  (S+SI)  and  helper  T-cells  (H+HI)  implied  by  the 
DeLisi  model  and  shows  the  absence  of  feedback  between  these  cells  required 
to  model  immune  responses.  For  this  reason  this  model  will  not  be  treated  in 
further  detail  in  this  paper. 

The  circuits  in  influence  diagrams  are  important  in  the  treatment  of  this 
paper.  A  circuit  in  an  influence  diagram  consists  of  a  path  which  starts  with 
a  given  vertex  and  follows  various  edges  in  the  directions  of  the  arrows  until 
the  original  vertex  is  reached  again.  The  length  of  a  circuit  is  the  number  of 
edges  in  the  circuit  that  must  be  traversed  from  a  given  vertex  until  the  same 
vertex  is  reached  again.  A  circuit  is  negative  if  it  has  an  odd  number  of  negative 
edges  and  positive  if  it  has  an  even  number  of  negative  edges  or  no  negative 
edges.  Following  previous  practice  (King,  1980,  1982,  1983)  positive  and  negative 
circuits  of  length  n  will  be  called  Bn  and  Cn  circuits,  respectively.  The  bipartite 
graphs  discussed  in  this  paper  can  have  circuits  only  of  even  lengths. 


Now  consider  the  dynamics  of  a  system  depicted  by  a  given  influence  diagram 
as  represented  by  the  flow  topology  around  an  unstable  equilibrium  point.  Possible 
flow  topologies  include  one  or  more  stable  steady  states  as  well  as  oscillations 
or  chaos.  The  earlier  work  on  chemical  reaction  networks  focused  on  networks 
exhibiting  periodic  oscillations  (King,  1980)  or  aperiodic  chaos  (King,  1983).  The 
current  work  on  immunological  control  networks,  on  the  other  hand,  is  concerned 
with  networks  proceeding  from  an  initial  virgin  state  by  exposure  to  antigen 
to  give  one  or  more  steady  states  corresponding,  for  example,  to  suppression 
or  help  configurations  (Herzenberg,  Black,  Herzenberg,  1980).  Suitable  influence 
diagrams  for  modelling  immunological  control  networks  thus  exhibit  an  appropriate 
number  of  stable  configurations  rather  than  oscillation  or  chaos. 

The  procedure  for  determining  the  flow  topology  corresponding  to  a  given 

influence  diagram  uses  methods  based  on  switching  circuit  theory  (Caldwell, 

1958;  Hu,  1968)  as  adapted  by  Glass  and  Kauffmann  (1973),  Glass  (1975a,  b,  1977a, 

b),  and  Glass  and  Pasternak  (1978)  and  which  have  been  termed  kinetic  logic 

> 

(Thomas,  1979).  In  the  logical  description  used  in  this  approach  time  may  be 
involved  either  synchronously  (Glass  and  Kauffmann,  1973;  Glass,  1975)  or 
asynchronously  (Thomas,  1984;  Kaufman,  Urbain,  and  Thomas,  1985).  Since  the 
special  features  of  the  more  complicated  asynchronous  version  (e.g.,  the  effects 
of  time  delays)  are  not  necessary  for  this  work,  the  simpler  synchronous  version 
will  be  used  throughout  this  paper.  For  simple  influence  diagrams  such  as  those 
discussed  in  this  paper  both  versions  appear  to  give  similar  flow  topologies. 

Consider  the  immunological  control  network  represented  by  an  influence 
diagram  as  a  synchronous  switching  network  in  which  time  is  quantized  so  that 
the  signs  of  the  first  time  derivatives  of  the  internal  variables  at  time  t+1  are 
determined  by  their  signs  at  time  t  (Glass,  1975a).  The  switching  state  at  any 
time  of  such  a  network  containing  n  internal  variables  can  be  represented  by 
an  n-dimensional  vector  of  Ts  and  0's  corresponding  to  positive  and  negative 


time  derivatives,  respectively,  of  the  n  internal  variables.  Such  an  n-vector 
is  called  a  state  vector.  The  number  of  different  possible  state  vectors  in  an 
n-dimensional  system  is  2n  and  these  state  vectors  can  be  represented  by  the 
2n  vertices  of  an  n-dimensional  cube.  The  possible  transitions  from  states  at 
synchronous  time  t  to  those  at  time  t+1  may  then  be  represented  by  arrows 
directed  along  the  edges  of  the  n-dimensional  cube.  In  such  transitions  the  value 
of  exactly  one  component  of  the  state  vector  changes.  Furthermore,  in  this 
synchronous  treatment  the  discrete  time  scale  t  is  chosen  so  that  it  advances 
one  unit  each  time  a  single  component  of  tho  state  vector  changes.  The  resulting 
n-dimensional  cube  with  directed  edges  is  called  a  state  transition  diagram.  The 
center  of  the  n-dimensional  cube  representing  a  state  transition  diagram 
corresponds  to  a  steady  state  in  which  all  of  the  first  time  derivatives  of  the 
internal  variables  are  zero  and  therefore  do  not  correspond  to  either  Boolean 
variable  0  or  1.  If  this  steady  state  is  unstable,  the  transitions  represented  by 
the  directed  edges  of  the  state  transition  diagram  define  the  fundamental  topology 
of  the  flow  in  the  neighborhood  of  the  unstable  steady  state. 

The  calculation  of  a  state  transition  diagram  corresponding  to  a  given  turbulent 
influence  diagram  can  be  performed  by  the  following  five-step  procedure: 

(1)  A  logical  relation  (OR  or  AND)  is  assigned  between  each  pair  of  edges  directed 
towards  a  given  turbulent  vertex.  An  OR  relationship  between  arrows  from 
vertices  X  and  Y  both  going  to  a  turbulent  vertex  Z  implies  that  the  internal 
variables  X  and  Y  affect  the  level  of  Z  in  separate  processes  corresponding  to 
an  equation  of  the  general  type 

dZ/dt  =  aXr  +  bYs  (1) 

An  AND  relationship  between  arrows  from  vertices  X  and  Y  both  going  to  a 


turbulent  vertex  Z  implies  that  the  internal  variables  X  and  Y  affect  the  level 
of  Z  in  the  same  process  corresponding  to  an  equation  of  the  type 

dZ/dt  =  aXrYs  (2) 

The  OR  relationship  at  turbulent  vertices  corresponding  to  separate  processes 
.appears  to  be  more  realistic  for  influence  diagrams  based  on  immunological  control 
networks  (Kaufman,  Urbain,  and  Thomas,  1985)  and  is  used  in  this  paper.  In  the 
case  of  influence  diagrams  modelling  chemical  reaction  networks  either  OR 
or  AND  relationships  at  turbulent  vertices  may  be  realistic  depending  upon  the 
nature  of  the  chemical  system. 

(2)  A  local  truth  table  is  generated  for  each  circuit  in  the  influence  diagram 
to  indicate  possible  transitions  between  state  vectors.  In  order  to  see  how  such 
a  local  truth  table  is  generated,  consider  possible  effects  of  one  internal  variable 
X  on  a  second  internal  variable  Y.  If  X  helps  Y  as  indicated  by  a  positive  arrow 
from  X  and  Y  in  the  circuit  under  consideration,  then  in  the  local  truth  table 
the  values  for  Y  in  each  possible  state  vector  at  time  t+1  will  correspond  to 
the  values  for  X  at  time  t.  If  X  suppresses  Y  as  indicated  by  a  negative  arrow 
from  X  to  Y  in  the  circuit  under  consideration,  then  in  the  local  truth  table  the 
values  fo  Y  in  each  possible  state  vector  at  time  t+1  will  be 'the  opposite  of 
the  values  for  X  at  time  t  (i.e.,  a  0  for  X  at  time  t  leads  to  a  1  for  Y  at  time 
t+1  and  vice  versa). 

(3)  The  local  truth  tables  for  each  of  the  circuits  in  the  influence  diagram  are 
combined  to  give  a  global  truth  table  using  the  logical  relationships  at  the  turbulent 
vertices  as  follows: 

(a)  If  any  of  the  local  truth  tables  show  a  value  of  0  for  a  variable  represented 
by  an  AND  turbulent  vertex,  then  the  global  truth  table  for  that  variable  also 


shows  a  value  of  0. 


(b)  If  any  of  the  local  truth  tables  show  a  value  of  1  for  a  variable  represented 
by  an  OR  turbulent  vertex,  then  the  global  truth  table  for  that  variable  also 
shows  a  value  of  1. 

(4)  Each  state  vector  in  the  global  truth  table  at  time  t+1  is  compared  with  the 
corresponding  state  vector  at  time  t.  If  any  such  pair  of  state  vectors  are  the 
same,  then  that  state  vector  corresponds  to  a  stable  configuration  of  the 
immunological  control  network  modelled  by  the  influence  diagram. 

(5)  If  no  stable  configurations  are  found,  then  possibilities  for  oscillatory  and/or 
chaotic  behavior  are  explored  by  using  the  global  truth  table  to  determine  the 
corresponding  state  transition  diagram  as  discussed  in  a  previous  paper  (King, 
1983).  Since  this  situation  does  not  arise  in  the  immunological  control  networks 
discussed  in  detail  in  this  paper,  the  relevant  algorithm  will  not  be  repeated 
here. 


3.  The  Herzenberg  Models 


The  Herzenberg  models  for  immunological  control  networks  (Herzenberg, 
Black  and  Herzenberg,  1980)  use  five  independent  variables:  B  cells,  two  types 
of  helper  T  cells  Thl  and  Th2,  and  two  types  of  suppressor  T  cells  Tsl  and  Ts2. 
Helper  Th  cells  carry  antiidiotypic  V|_|  structures  and  thus  correspond  to  negative 
vertices  in  the  bipartite  immunological  control  networks  (ellipses  in  Figure  2). 
Suppressor  Ts  cells  and  B  cells  carry  idiotypic  V|_ )  structures  which  are 
complementary  to  the  antiidiotypic  Vf_j  structures  and  thus  correspond  to  positive 
vertices  in  the  bipartite  immunological  control  networks  (rectangles  and  squares 
in  Figure  2).  The  soluble  products  of  the  Herzenberg  models  (e.g.,  P-\,  P2,  and 
P3  in  Eisenfeld  and  DeLisi,  1985)  are  viewed  as  dependent  rather  than  independent 


variables  with  their  levels  being  so  tightly  coupled  to  B  cells  that  they  do  not 
need  to  be  considered  as  separate  variables.  In  this  way  the  number  of  variables 
in  the  Herzenberg  models  considered  in  this  paper  can  be  limited  to  five. 

The  Herzenberg  models  consist  of  a  bistable  core  regulatory  circuit  (CRC) 
which  can  be  maintained  stably  in  either  the  help  or  suppression  mode.  One 
or  more  auxiliary  regulatory  circuits  (ARC's)  are  linked  to  this  CRC  and  can 
switch  the  CRC  between  the  help  and  suppression  modes. 

The  Herzenberg  CRC  is  depicted  in  Figure  2  (top).  It  represents  the  simplest 
way  of  linking  helper  and  suppressor  T  cells  to  satisfy  the  following  conditions: 

(1)  Antiidiotypic  cells  can  only  interact  with  idiotypic  cells  or  vice  versa  so  that 
only  bipartite  graphs  are  allowed  thereby  excluding  graphs  containing  circuits 
of  odd  lengths. 

(2)  The  Th  and  Ts  cells  are  differentiated  so  that  the  populations  of  target  Th 
and  helper  Th  for  a  given  Ts  are  drawn  from  different  Th  populations. 

The  Herzenberg  CRC  (Figure  2,  top)  is  thus  seen  to  consist  of  a  single  positive 
circuit  of  length  4,  called  a  fundamental  quadrilateral  (King,  1982).  By  a  theorem 
in  the  previous  paper  (Theorem  8,  in  King,  1982),  a  positive  fundamental 
quadrilateral  must  have  two  attracting  regions  or  stable  configurations  exactly 
as  suggested  in  the  Herzenberg  analysis. 

In  order  to  determine  the  nature  of  these  two  stable  configurations  the  truth 
table  for  the  Herzenberg  CRC  is  determined  (Table  1).  The  stable  "help" 
configuration  corresponds  to  a  (Tsl,  Thl,  Ts2,  Th2)  state  vector  of  (0,  1,  1,  0) 
or  dTsl/dt  <  0,  dThl/dt  >  0,  dTs2/dt  >  0,  and  dTh2/dt  <  0.  In  this  configuration 
Thl  increases  (dThl/dt  >  0)  and  thus  helps  the  differentiation  of  Ts2  (i.e., 
dTs2/dt  >  0)  which  then  depletes  Th2  (i.e.,  dTh2/dt  <  0),  which  in  turn  disables 
differentiation  and  expression  of  the  Ts2  population  (i.e.,  dTs2/dt  <  0),  which, 
if  present  and  active,  would  be  capable  of  attacking  Thl.  The  stable  "suppression" 


configuration  corresponds  to  a  (Tsl,  Thl,  Ts2,  Th2)  state  vector  of  (1,  0,  0,  1) 
or  dTsl/dt  >0,  dThl/dt  <  0,  dTs2/dt  <  0,  and  dTh2/dt  >  0  leading  analogously 
to  an  increase  in  Th2  and  a  decrease  in  Thl. 

The  shifting  of  the  CRC  between  the  help  and  suppression  configurations 
can  be  accomplished  by  addition  of  an  ARC.  The  simplest  ARC's  consist  of  a 
single  strong  vertex,  corresponding  to  the  B  cells  and  their  interactions  with 
the  Th  cells  through  soluble  products.  The  single  strong  ARC  vertex  can  be  linked 
to  the  CRC  in  two  ways  to  preserve  the  bipartite  nature  of  the  entire  network 
(Figure  2,  middle  and  bottom).  The  first  of  these  networks,  B4  +  C4,  (Figure 
2,  middle)  corresponds  to  the  "suppressive  ARC"  Herzenberg  model  (Herzenberg, 
Black  and  Herzenberg,  1980:  Scheme  2,  p.5)  and  consists  of  two  circuits  of  lengths 
4  of  opposite  signs  with  two  edges  in  common  and  a  single  turbulent  vertex,  namely 
Th2.  This  turbulent  vertex  is  assigned  the  OR  relationship  since  the  interactions 
of  B  and  Ts2  with  Th2,  as  indicated  by  the  directed  edges,  are  separate  processes. 
The  truth  table  for  this  B4  +  C4  network  (Table  2)  illustrates  the  process  of  first 
generating  local  truth  tables  for  each  circuit  and  then  combining  the  local  truth 
tables  into  a  global  truth  table  by  using  the  logical  relationships  at  turbulent 
vertices  (King,  1983).  From  the  global  truth  table  in  Table  2  only  one  of  the 
32  possible  (Tsl,  Thl,  Ts2,  Th2,  B)  state  vectors,  namely  (1,  0,  0,  1,  0),  is  seen 
to  lead  to  a  stable  configuration.  This  stable  configuration  has  dTsl/dt  >  0, 
dThl/dt  <  0,  dTs2/dt  <  0,  and  dTh2/dt  >  0  corresponding,  to  the  stable  "suppression" 
configuration  of  the  isolated  CRC  (Table  1).  This  simple  example,  which 
corresponds  to  one  of  the  Herzenberg  networks,  shows  how  addition  of  an  ARC 
to  a  CRC  can  establish  the  CRC  into  one  of  its  two  stable  configurations. 

The  other  possible  way  of  adding  a  single  strong  vertex  (B)  to  the  Herzenberg 
CRC  while  preserving  the  bipartite  nature  of  the  network  (Figure  2,  bottom) 
leads  to  a  five-vertex  B4  +  B2  network  consisting  of  positive  circuits  of  length 


4  and  2  sharing  a  common  vertex.  This  shared  vertex  (Thl)  is  thus  turbulent 
and  is  assigned  the  OR  relationship  as  in  the  64  +  C4  case  above.  The  truth  table 
for  this  84  +  B2  network  (Table  3)  indicates  a  single  stable  configuration 
corresponding  to  the  stable  "help"  configuration  of  the  isolated  CRC  (Table  1) 
with  dTsl/dt  <  0,  dThl/dt  >  0,  dTs2/dt  >  0,  and  dTh2/dt  <  0  opposite  to  the  stable 
"suppression"  configuration  of  the  B4  +  C4  network  discussed  above. 

This  analysis  thus  indicates  simple  ways  of  attaching  a  fifth  vertex  (B)  to 
a  given  vertex  (Thl)  of  the  simple  four-vertex  Herzenberg  CRC  to  lock  it  in 
either  the  "suppression"  or  "help"  configurations.  The  help-stimulating  Herzenberg 
network  (Herzenberg,  Black  and  Herzenberg,  1980:  Scheme  3,  p.  6)  is  a  sum  of 
the  B4  +  C4  and  B4  +  B2  networks  in  Figure  2  (middle  and  bottom)  with  the  B4 
+  B2  behavior  obviously  dominating  (i.e.,  a  greater  weight  of  the  B-*-Th1  edge 
relative  to  the  B-*Th2  edge).  The  allotype-suppressive  Herzenberg  network 
(Herzenberg,  Black  and  Herzenberg,  1980:  Scheme  4,  p.8)  is  still  more  complicated 
consisting  of  a  seven-vertex  B4  +  C4  +  C4  +  82  network  having  a  triply  turbulent 
Thl  vertex  common  to  all  four  feedback  circuits  as  well  as  a  (doubly)  turbulent 
Th2  vertex. 

4.  The  Model  of  Kaufman,  Urbain,  and  Thomas 

The  model  of  Kaufman,  Urbain  and  Thomas  (1985)  for  an  immunological  control 
network  uses  five  independent  variables:  B  =  B  cells,  H  =  Th  =  helper  T  cells, 

5  =  Ts  =  suppressor  T  cells,  A  =  antibody,  and  E  =  antigen  (epitope).  The 
relationship  between  these  five  variables  is  depicted  in  the  network  at  the  top 
of  Figure  3.  In  this  network  A  and  E  are  weak  vertices  and  B,  H,  and  S  are  strong 
vertices.  Removal  of  the  vertices  A  and  E  and  the  associated  five  edges  E  -*•  B, 
E -►  H,  E  *A,  B  *  A,  and  H  -*■  A  leads  to  a  network  with  the  three  vertices  B, 


H,  and  S.  In  this  stage  B  becomes  a  weak  vertex  because  of  the  removal  of  the 
only  edge  directed  to  B,  namely  E-*-  B.  Thus  a  second  stage  involving  the  removal 
of  the  vertex  B  and  associated  edge  B  -*•  H  is  necessary  to  obtain  a  two-vertex 
influence  diagram  corresponding  to  the  original  five-vertex  immunological  control 
network. 

At  this  point  there  arises  the  question  of  treating  the  loops  at  the  vertices 
H  and  S  representing  self-interactions.  Deletion  of  these  loops  leading  to  a  simple 
C2  circuit  (Figure  3,  bottom  left)  is  the  simplest  approach  but  results  in  the 
loss  of  important  information.  The  flow  topology  of  such  a  simple  C2  circuit 
corresponds  to  limit  cycle  oscillation  but  such  a  system  can  also  function  as 
a  biochemical  switch  (Glass  and  Kauffman,  1973).  A  more  complicated  but  better 
approach  is  the  replacement  of  these  self-interaction  loops  with  a  positive  circuit 
between  two  complementary  idiotypes  (Figure  3,  bottom  right).  This  influence 
diagram  has  the  same  four  variables  as  the  Herzenberg  CRC  (Figure  2,  top). 
However,  these  variables  are  not  only  linked  in  different  ways  but  also  must 
be  assigned  different  parities  to  preserve  the  bipartite  nature  of  this  graph.  The 
feedback  circuit  structure  of  the  influence  diagram  arising  from  the  model  of 
Kaufman,  Urbain  and  Thomas  (1985)  is  much  more  complicated  than  the  simple 
84  circuit  of  the  Herzenberg  CRC.  As  a  result  of  this,  all  four  vertices  of  the 
influence  diagram  (Table  3,  bottom  right)  arising  from  the  model  of  Kaufman, 
Urbain  and  Thomas  (1985)  are  turbulent. 

The  truth  table  (Table  4)  corresponding  to  this  influence  diagram  (Figure 
3,  bottom  right)  can  be  determined  by  assigning  OR  relationships  to  each  of  the 
turbulent  vertices.  This  leads  to  two  stable  configurations  as  was  found  for  the 
analogous  analysis  of  the  Herzenberg  CRC  discussed  above.  One  of  the  stable 
configurations  corresponds  to  a  (HI,  H2,  SI,  S2)  state  vector  of  (0,  0,  1,  1)  or 
dH1/dt  <  0,  dH2/dt  <  0,  dSI/dt  >  0,  and  dS2/dt  >  0.  This  corresponds  to  the 
"non-responsive"  state  of  Kaufman,  Urbain,  and  Thomas  (1985)  in  which  only 
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Ts  cells  are  present.  The  other  stable  configuration  corresponds  to  a  (HI,  H2, 
SI,  S2)  state  vector  of  (1,1, 1,1)  or  dH1/dt  >  0,  dH2/dt  >  0,  dSI/dt  >  0,  and  dS2/dt 
>  0.  This  corresponds  to  the  "memory"  state  of  Kaufman,  Urbain,  and  Thomas 
(1985)  in  which  both  Th  and  Ts  cells  are  present.  Thus  this  treatment  of  the 
Kaufman,  Urbain  and  Thomas  model  agrees  with  the  results  of  their  treatment 
for  the  case  of  the  absence  of  antigen  (E).  Since  the  weak  antigen  vertex  (E) 
was  deleted  from  the  original  network  (Figure  3,  top)  in  deriving  the  corresponding 
influence  diagram  (Figure  3,  bottom  right),  this  is  the  case  to  which  the  calculation 
in  Table  4  corresponds.  This  also  demonstrates  the  essential  equivalence  of  the 
method  in  this  paper  to  the  method  used  by  Kaufman,  Urbain,  and  Thomas  (1985) 
in  a  case  where  either  method  can  be  applied. 

5.  Other  Models 

Two  older  models  for  immunological  control  networks  have  also  been  examined 
by  the  methods  outlined  in  this  paper.  The  model  of  Hiarnaux  (1977)  consists 
of  a  cycle  of  C2  circuits  of  the  antibodies  Abk  (1  _<  k  <_  n)  with  an  "external"  C2 
circuit  to  the  antigen  (Ag)  (Figure  4).  Investigation  of  the  simplest  non-trivial 
bipartite  case  where  n  =  4  and  using  OR  relationships  at  the  four  turbulent  vertices 
Abk  (1  ±  k  <_  4)  indicates  only  one  stable  configuration  corresponding  to  an  (Ag, 
Abl,  Ab2,  Ab3,  Ab4)  state  vector  of  (0,  1,  1,  1,  1).  The  model  of  Richter  (1978a, 
b)  consists  of  an  open  chain  of  C2  circuits  (Figure  5).  Investigation  of  the  flow 
topology  of  this  system  for  n  *  4  reveals  the  absence  of  any  stable  configurations 
suggesting  oscillatory  behavior  of  this  model.  Hoffmann  (1982)  has  already 
indicated  difficulties  of  these  models  relating  to  the  absence  of  the  required 
stable  configurations. 
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6.  Summary 

This  paper  shows  that  flow  topological  methods  developed  for  the  analysis 
of  the  dynamics  of  chemical  reaction  networks  (King,  1980,  1982,  1983;  Class 
and  Kauffmann,  1973;  Glass,  1975a,  b,  1977a,  b;  Class  and  Pasternak,  1978)  can 
also  be  used  for  the  analysis  of  immunological  control  networks  with  appropriate 
modifications  to  account  for  the  differences  between  the  two  types  of  networks. 
Thus  the  chemical  reaction  networks  of  interest  have  included  systems  exhibiting 
periodic  oscillations  or  aperiodic  chaos  whereas  the  immunological  control 
networks  represent  nonoscillatory  systems  having  one  or  more  stable  configurations 
(e.g.,  "help"  and/or  "suppression").  In  addition,  the  complementary 
idiotype-antiidiotype  interactions  in  the  immunological  control  networks  lead 
to  bipartite  graphs,  a  feature  not  generally  found  in  chemical  reaction  networks. 

The  Herzenberg  models  for  immunological  control  networks  (Herzenberg, 
Black,  and  Herzenberg,  1980)  are  most  suitable  for  the  methods  outlined  in  this 
paper  because  of  the  presence  of  a  minimum  number  of  turbulent  vertices.  Flow 
topological  analysis  of  the  Herzenberg  models  leads  to  the  locking  in  stable  "help" 
and/or  "suppression"  configurations  observed  experimentally.  In  addition,  the 
flow  topology  of  the  model  of  Kaufman,  Urbain,  and  Thomas  (1985)  confirms 
the  presence  of  the  "memory"  and  "non-responsive"  states  found  in  their  analysis 
indicating  essential  agreement  of  the  two  methods  for  the  treatment  of  cases 
where  both  can  be  applied.  However,  the  flow  topologies  of  the  models  of  Hiernaux 
(1977)  and  Richter  (1978a,  b)  indicate  the  presence  of  only  one  and  zero  stable 
configurations,  respectively,  thereby  suggesting  major  difficulties  with  these 
models  in  accord  with  their  recent  analysis  of  Hoffmann  (1982).  In  any  case 
this  paper  demonstrates  the  application  of  flow  topology  for  the  analysis  of 
immunological  control  networks  to  reproduce  some  of  the  basic  features  of  such 


networks  in  simple  systems  using  a  few  as  four  or  five  internal  independent 
variables. 
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TRUTH  TABLE  FOR  THE  MODEL  OF  KAUFMAN,  URBAIN,  AND  THOMAS 


Figure  1  :  The  immunological  control  network  of  DeLisi  (1983).  (a)  Top:  the 
original  seven-vertex  directed  graph;  (b)  Middle:  removal  of  the  weak 
vertices  A  and  B  and  the  associated  edges  A  -*■  SI,  A  -*•  HI,  H  -*■  B,  and 
M$-*-  B;  (c)  Bottom:  the  corresponding  two-vertex  signed  directed  graph 
in  which  the  (weak)  vertices  S  +  SI  and  H  +  HI  correspond  to  suppressor 
and  helper  T  cells,  respectively. 

Figure  2:  Influence  diagrams  corresponding  to  the  Herzenberg  models  (Herzenberg, 
Black  and  Herzenberg,  1980).  (a)  Top:  the  four-vertex  CRC;  (b)  Middle: 
addition  of  a  fifth  vertex  (B)  to  give  the  B4  +  C4  system  stabilizing 
the  "suppression"  configuration;  (c)  Bottom:  addition  of  a  fifth  vertex 
(B)  to  give  the  B4  +  82  system  stablizing  the  "help"  configuration. 

Figure  3:  The  immunological  control  network  of  Kaufman,  Urbain,  and  Thomas 
(1985).  (a)  Top:  the  original  five-vertex  signed  directed  graph;  (b)  Upper 
middle:  removal  of  the  weak  vertices  E  and  A  as  well  as  the  loops  at 
H  and  S  to  give  a  three-vertex  signed  directed  graph;  (c)  Lower  middle: 
removal  of  the  weak  vertex  B  and  the  edge  B  -*•  H  from  the  three-vertex 
signed  directed  graph  to  give  a  two-vertex  signed  directed  graph;  (d) 
Bottom  left:  removal  of  the  positive  loops  at  H  and  S;  (e)  Bottom  right: 
replacement  of  the  positive  loops  at  H  and  S  with  positive  circuits 
between  the  two  complementary  idiotype  pairs  H1/H2  and  S1/S2 
respectively. 
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Figure  4:  The  immunological  control  network  of  Hiernaux  (1977):  fa)  Top:  the 
general  case;  (b)  Bottom:  the  simplest  non-trivial  bipartite  case  of  this 
network  where  n  =  4  with  OR  relationships  at  the  four  turbulent  vertices 
Abk  (1  <_  k  <  4). 

Figure  5:  The  immunological  control  network  of  Richter  (1978). 
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